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ABSTRACT 


Two subjects are discussed, which are believed relevant to the structural 
analysis of vertical-axis wind turbines. The first involves the derivation 
of dynamic differential equations, suitable for studying the vibrations of 
rotating, curved, slender structures. The Hamiltonian procedure is advocated 
for this purpose. Various reductions of the full system are displayed, which 
govern the vibrating Troposkien when various order-of-magnitude restrictions 
are placed on important parameters. 

The final section discusses the possible advantages of the WKB asymptotic 
method for solving these classes of problems. A special case of this method 
is used illustratively to calculate eigenvalues and eigenfunctions for a 
"flat” turbine blade with small flexural stiffness. 

PRINCIPAL SYMBOLS 


k 


s 

t 



X = X +5 
o 

y = y^+ n 


Z = ? 


Argument of elliptic integrals & functions 
Length of blade 

(Constant) mass per unit length of blade 
Coordinate measuring distance along blade 
Time 

Limits in Hamiltonian integral 


Instantaneous coordinates of point on vibrating blade 


Xo(s), y^(s) 
X 

m 


El 


s 


Coordinates of blade in rotational equilibrium 
Coordinate of blade end 
Maximum value of y 

o 

Flexural rigidity for bending in x-y-plane 
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K, K 

o 

L 

T= T + T 

U(s) 

6 (• • • ) 

X 

y 

(JL) 

a 


Flexural rigidity for bending normal to x-y-plane 

(Constant) x-component of equilibrium tensile force 

Curvatures defined by Eqs. (4) 

Kinetic minus potential energy 
Instantaneous tensile force in blade 

Equilibrium value of T 

Vector of unknowns in WKB method 
Variation of a quantity 
Eigenvalue in WKB method 
Parameter defined by Eq. (9) 

Circular frequency of simple harmonic motion 
Angular velocity about support axis of blade 


Superscripts, Etc . 

( ) Partial derivative with respect to time. 

( ) ' Total derivative with respect to s 


INTRODUCTION 


The title of this paper is somewhat misleading, since it actually addresses 
two subjects relevant to the structural analysis of vertical-axis wind turbines. 
The first involves the derivation of differential equations needed for the 
vibration and aeroelastic analysis of a rotating, curved, slender structure 
like a beam in the Troposkien shape. It is emphasized that a systematic 
procedure, such as that furnished by Hamilton's principle, can be useful both 
for ensuring that the results are correct and for developing a hierarchy of 
simplifications when various approximations are made. 

The second subject concerns the possible advantages of what Steele (Ref. 1) 
chooses to call the WKB (Wentzel-Kramers-Brillouin) method for solving certain 
of these problems. An elementary Illustrative example is offered. 

Admittedly, detailed structural design studies of Darrieus-type wind 
turbines are best carried out by means of finite-element approximations. This 
approach is unavoidable when the designer is faced with such complications 
as interactions with the support structure, reinforcing struts, or variable 
properties along the blades. Calculations of operating stresses, vibration 
modes, etc. by finite elements appear in several published papers. Examples 
are Weingarten and Nickell (Ref. 2); Weingarten and Lobitz (Ref. 3); and 
Biffle (Ref. 4). Many other related citations are given in the survey paper 
by Blackwell et al . (Ref. 5). 

Parametric studies and the goal of fundamental understanding are not so 
well served, on the other hand, by the purely numerical investigations of point 
designs. It is believed that Troposkien vibrations offer some interesting 
problems in applied mechanics, and the discovery of analytical or semi-analytical 
solutions to appropriate differential equations can certainly not be ruled out. 
The aim of this paper is to encourage that search. 
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Figure 1 depicts a Darrieus machine, with a typical blade both in its 
unperturbed equilibrium shape of YqCxq) and slightly deformed. The instan- 
taneous time-dependent position of mass element xK ds is specified by coordinates 
x,y,z, measured in a frame rotating with the constant angular velocity . 
Various simplifications are here adopted from the outset. Thus, attention is 
focused on a single, uniform blade by assuming its ends to be fully restrained 
at the points (i Xjj^, o, o). Gravitational force is neglected ( cf . comparative 
studies cited by Weingarten and Lobitz Ref. 3). 


The current experimental Darrieus configurations deviate somewhat from 
the true Troposkien shape (Blackwell and Reis, Ref. 6), and it is known that 
these deviations can produce substantial equilibrium bending stresses (Ref. 3). 
Here yo(xo) is chosen, however, to be the perfect "skiprope." This is a 
restriction that is easy to remove, as are the prescribed limitations to 
essentially infinite torsional and elongational rigidity of the blade. In fact, 
during free vibration, torsion is often nearly uncoupled from the other degrees 
of freedom when the blade elastic axis and line of C.G.’s coincide — a design 
feature that would seem desirable in practice. 

VIBRATION EQUATIONS FOR SEVERAL IDEALIZED TROPOSKIENs" ^ 


It is well known (e.g.. Ref. 6) that the zero bending moment shape is 
governed by 






( 1 ) 


With support conditions as in Fig. 1, the solution of (1) is conveniently 
expressed as 



sn 


K(k) I 1 + ^ j ; k 


( 2 ) 


Here sn is Jacobi’s elliptic function, whereas K is the complete elliptic 
integral of the second kind with argument 


k = 


1 + 


4H 

o 

mU y 
•'m 


-is 


__ 2 2 

Equation (3) introduces a dimensionless group ^ ^ which 

slight modifications, is perhaps the key parameter for vibration 
studies. 


(3) 

, with various 
and stability 


* 

Principal symbols are defined in Fig. 1 and the list at the beginning. 
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As suggested by Mallett in his work on the catenary (Ref. 7), distance s 
along the blade seems to be an efficient choice for the spatial independent 
variable. In terms of s, the instantaneous curvatures of the undeformed and 
deformed shapes can be written, respectively. 



(4a) 


(4b) 


For the inextensible blade with flexural stiffnesses El and El , an appro- 

s ^ 


priate Hamiltonian is 
•^2 




-It 

1 


L + (Length constraint) dt = % 


] 


• 2 2 2 2 
+ 2Qyz - 2Qzy + Q y + z 


J I + 


(5) 


+ T 


(s>t) 1^1 j \9s) \9s/ J j 


dsdt 


It is easily seen that the Lagrange multiplier on the constant length constraint 
is simply the tension T(s,t). 

The principle 6H = 0 is enforced by taking variations of (5) on x,y,z 
and T. Suitable integrations by parts with respect to s or t then lead to the 
general nonlinear dynamical system: 



(6a) 


(6b) 


(6c) 


(6d) 
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Small perturbations n , C and x are introduced on the dependent 

variables, products of small quantities are neglected, and the equilibrium 
relations corresponding to (1) are subtracted out. There results the most 
general system, governing the free dynamics of a uniform, inextensible 
Troposkien; 



In principle, the assumption of simple harmonic time dependence at frequency a> 
converts Eqs. (7) into an ordinary system, whose eigenvalues and eigenfunctions 
describe the free vibration. It is also not difficult to add forcing terms, 
such as aerodynamic loads. 

Equations (7), as written, offer little attraction to the analyst. Accord- 
ingly, several simpler versions will be discussed, along with the approximations 
that lead to them. For reasons of space, little attention is paid to boundary 
conditions, but their formulation is not a difficult matter. 

!• The Rotating Rope or Chain 

When El = El = 0, a catenary-like structure is left. Its linearized 

s z 

differential equations are 


m 5 

3 

3s 


II 

0 

m l^n 

- 2 


In 

9s 


m 

+ 2nfi 


T 

o 

9? 1 
3s J 


(8a) 

= 0 (8b) 

0 , (8c) 


with (7d) unchanged. Equations (8) are not relevant to wind turbines. The 
reason for displaying them is to make the point that, even in this case, the 
four variables remain coupled. In addition to the familiar centrifugal term 
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> Coriolis effects couple the in- and out-of-plane degrees of free om. 

These latter gain in importance relative to the centrifugal as the vi ra ion 
frequency increases. 

When (8) are made dimensionless, a two-parameter system is found, involving 

_ 2 2 

_ mo) (9) 

y = -H 

o 

2 2 

and the frequency ratio fi/u . It appears that y /w is 

interesting designs, fi/u is less than unity, but it remains ( ) or e 

lower vibration modes. 

2. Small In-Plane and Large Out-of-Plane Bendin g Stiffnesses 

It is obvious that Darrieus turbine blades with typical airfoil shapes will 

have 


El 
z 

El 

s 


» 1 


( 10 ) 


Indeed, Ref. 3 quotes a value over 40 for this ratio on the SANDIA 17 meter 
design. An interesting reduction of system (7) is found when both inequality 
(10) is satisfied and when 


El 


H X 

o m 


« 1 


( 11 ) 


(this ratio is estimated to be 0.015 for the 17-meter). It can then be 
reasoned that ? « n , so that the Coriolis term in (7b) is negligible, and 
that the flexural terms may be dropped from (7a, b). The resulting equations are 


m ? - 

3s 

["o 

3s 

+ Xq ’ =0 

(12a) 

m l^n - 



T 

o 

+ y ' t1 = 0 

3s ^ o J 

(12b) 

3s 

dx 

3n _ 
3s 

0 


(12c) 


o 
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Three coupled dependent variables are involved in (12). Nevertheless, 
they constitute a considerable simplification and might form the basis for 
useful parametric studies. At considerable labor, one can eliminate C(s,t) 
among (12 a,b,c). The result does not seem of great value, however, and 
efforts to achieve a single equation in one unknown where unsuccessful. 

3. The ”Flat" Blade with Small Bending Stiffness 

If in addition to the other assumptions behind system (12) the approximation 
is made that 


dx 


« 1 


(13) 


then dropping the t term in (12b) is justified. A further replacement of 
where is its large, constant x-component leads to 



3x 


2 


o 


m ds 

H dx 
o o 




0 


(14) 


This is essentially the in-plane-deformatlon equation employed by Ham (Ref. 8), 
in his investigations on free vibration and flutter of the rotating blade. 

In view of the apparent limitations on its validity, further study would seem 
justified into the parameter ranges within which it can be used in practice. 

DISCUSSION OF ASYMPTOTIC METHODS AND AN APPLICATION 

The excellent review by Steele (Ref. 1) removes the need for any recapit- 
ulation here of how these methods are applied to ordinary differential equations 
arising in solid mechanics. The particular form useful in vibration problems 
with spatially-varying coefficients is known by such names as "phase integral," 
WKB and WKBJ. In the form most broadly applicable here, it starts with a 
homogeneous problem formulation in terms of a "state vector" U(s): 


= A(s;X) U(s) (15) 

Here A(s,X) is a prescribed coefficient matrix, containing an eigenvalue X 
which, in some sense, is a large quantity. A is expanded according to 

A (s;X) = XA^(s) + Aj^(s) + ^. 2 ( 5 ) + ••• (16) 

The series in (16) may be convergent or asymptotic. In either event, clever 
changes of variable are employed to construct a corresponding succession of 
approximate solutions. 


45 



After the assumption of simple harmonic motion, it is worth noting that all 
the systems presented in the preceding section can be recast in state vector 
form. The same is true of equations for flutter stability or forced sinusoidal 
response. The principal recommendation of this paper is that the WKB theory 
should be put to work in order to determine what contribution it may make to 
the subject at hand. 


All that can be offered at present is a simple example to illustrate the 
possibilities. For harmonic oscillations at frequency o) , with mode shape 

tT(x ), it is an easy matter to recast the ’’flat - Troposkien" eq. (14) as follows 
o 


dx^ 


Ail r 

L 


1 + 


n2 2 
mU y 
•^m 

H 


1 - 




m 


n = 0 , 


(17) 


where 


A 


+ 


(18) 


In view of (2), the quantity in brackets of (17) may be rewritten 

2 


[•••]■ 


- o2 2 

m y 

1 . 

1 + — n cn 


H 


(km (l + ^ r) 


(19) 


.Now (17), with the associated boundary conditions n (±x^) - 0 , is a 

special case of the problem solved by Steele (Ref. 1, Sect. 2), using the 
Green-Liouville transformation. In more general terms, Steele transforms 
the differential equation 




( 20 ) 


by introducing 



4 = ?(x) , 

y = \()(x) n(4) 

(21a,b) 

With 

the following 

optimal choice of these functions: 



(4’)^ = p/q 


(21c) 

and 

'P = ^pqj , 


(21d) 

(20) 

is reduced to 





a ] " = ° ’ 

(22) 


ds^ I- 

q J 
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where 


v(x) = r - 


(23) 


When X is large, the limiting form of the solution can be re-transformed to 


y(x) 


where 


C(x) 



(24) 


(25) 


has the name "phase integral." x, is a point associated with one of the 
boundary conditions. Eigenvalues are readily estimated by combining (24) 
(25) and the second B.C. 


and 


In the present example, one associates Steele's x variable with x /x 
y with n . Here also ° ™ 


P(x) 
q(x) 
r (x) 


= 1 


= 1 + 


m 



H 
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= 0 

-,2 

m A X 
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2 


H 

o 



(26a) 

(26b) 

(26c) 

(26d) 


After specialization of the general solution, it is found that a first approxi 
mation to the f ree-vibration eigenvalues is 

^n " > (n = 1,2, •••) , (27a) 

L 


where [cf. (26b)] 



The mode shape corresponding to X is 



(27b) 


(28a) 
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where 



(28b) 


It is even possible, by one additional integration, to obtain a second 
to th; eigpnvalupp T in ancordanca with > 
in Sect. 3.3, Ref. 1. Details are not given here, since correcti 
generally less than 5%. 

In view of (18), (26d) and other knowledge of the p 

Darrleus turbines, one can conclude that » is quite a large 

possibly for the fundamental solution of (17). According y, n 9 m n 8 

Ts Spued for values of the parameter k [eq. (3) ranging from 0 2 to 0.8. 
fk = 0 57 characterizes a turbine with diameter equal to its heig •; 

Table I Lrnishes some numerical results* for the first five eigenvalues and 
the corresponding dimensionless natural frequencies. 


Figure 2 plots some representative vibration mode shapes for k = 0.2. 
This parameter was chosen because the condition (13), implied in (17), is most 
likely to be satisfied accurately in this rather "flat" case. The mode 
corresponding to n= 1 is not shown, because it involves a low-frequency, 
symmetrical stretching motion. Although not ruled out by (17), it obviously 
violates the condition ( 6 d) of inextensibility and is regarded as a physical 
impossibility. For the same reason, n= 1 eigenvalues are enclosed by 
parentheses in Table I. 


Time has not permitted extensive comparisons between the foregoing results 
and those of previous finite-element analyses. One possibility is furnished, 
however, by the data in Fig. 5 of Ref. 3. These relate to a blade like those 
of the SANDIA 17-meter turbine, but with the stiffening struts removed; the 
rotation rate is 75 RPM. From the modal symmetries, it is clear that their 
numbers are equivalent to ( 0 ) 2 / 0 ) = 1.55 and ( 0 ) 3 / 0 ) - 2.95 at k - 0.5 . 
From the data in Table I, the corresponding numbers by the WKB approximation 
are 1.75 and 2.97, respectively. 


The present computations were carried out on the digital computer only 
for convenience. Nothing more is Involved than subroutines giving elliptic 
functions and the numerical evaluation of well-behaved Integrals. It is ^ 
believed unlikely that any much more efficient scheme exists for dealing with 
equations like (17). Further Troposkien dynamics investigations by means 
of the WKB method will hopefully be stimulated by this very elementary first 
attempt . 


*The author is indepted to Messrs. James Nathman and Larry Lehman for 
carrying out these calculations. 
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DISCUSSION 


Question: 

Since torsion is structurally coupled with out-of-plane bending in a curved 
beam, isn’t it true that whenever such bending occurs, torsion must occur? 

Comment 1 : Since coriolis forces dominate the coupling between in-plane and 

out-of-plane bending, blade chordwise elastic axis and center of gravity loca- 
tions need not coincide, and blade design effort and manufacturing cost are 
greatly reduced. 

Comment 2; In certain experimental cases of Darrieus blade flutter, in-plane 
and out-of-plane bending displacement perturbations are the same order of mag- 
nitude. 
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Answer; Agree with Comment 2 . Feel we must gain more experience with various 
sizes of rotor and other structural details before Comment 1 can be accepted 
as a design principle for all Darrieus machines. Helicopter rotors also have 
important Coriolis couplings, yet practice with them is to massbalance so that 
the C.G. and quarter-chord axes coincide. 

With regard to the question, author believes that the theory of his paper 
requires, as a strict condition, that the torsional rigidity GJ just be 
extremely large. He has found other possible circumstances, however, in whicn 
coupling with torsion may be negligible. This occurs when the C.G. is at the 
quarter chord and when a dimensionless parameter (involving GJ divided by tne 
mass moment of inertia in torsion) is large compare to unity. 


TABLE I. - SECOND-APPROXIMATE EIGENVALUES AND DIMENSIONLESS NATURAL 
FREQUENCIES FOR TROPOSKIEN BLADES OF SEVEN 

DIFFERENT ASPECT RATIOS 



— 

k 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

n 

y /x 
m m 

0.2625 

0.4100 

0.5807 

0.7908 

1.071 

1.486 

2.226 

1 

tl) 

X 

(1.511) 

(1.525) 

(1.489) 

(1.440) 

(1.376) 

(1.057) 

(1.181) 


(0.7395) 

(0.719) 

(0.689) 

(0.648) 

(0.593) 

(0.517) 

(0.406) 

2 

3.061 

2.957 

2.805 

2.599 

2.325 

1.964 

1.479 


z 

l.lhl 

2.168 

2.058 

1.904 

1.691 

1.388 

0.909 

3 

z 

% 

4.593 

4.444 

4.235 

3.966 

3.638 

3.251 

2,801 


Q) Q ^ 

3.545 

3.446 

3.306 

3.126 

2.907 

2.649 

2.356 

4 

< 

6. 124 

5.924 

5.642 

5.271 

4.797 

4.187 

3.360 


0), /n 

4.808 

4.677 

4.492 

4.247 

3.928 

3.507 

2.903 

5 

4 

% 

7.655 

7.405 

7.053 

6.593 

6.016 

5.304 

4.433 


U)^/ 

6.057 

5.895 

5.665 

5.364 

4.983 

4.510 

3.926 
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Figure 1. - Equilibrium and perturbed positions of a Troposkien shape 
which is vibrating while rotating about x^-axis with constant 
angular velocity si. 
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Figure 2. - First three physically possible natural modes of vibration 
for the Troposkien with (mx2ft2/H^)' = 1 .555, corresponding to k = 0.2. 
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